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Abstract 

We investigate the variation of the gas and the radiation pressure in accretion 
disks during the infall of matter to the black hole and its effect to the flow. While 
the flow far away from the black hole might be non-relativistic, in the vicinity of 
the black hole it is expected to be relativistic behaving more like radiation. There- 
fore, the ratio of gas pressure to total pressure (/3) and the underlying polytropic 
index (7) should not be constant throughout the flow. We obtain that accretion 
flows exhibit significant variation of f3 and then 7, which affects solutions described 
in the standard literature based on constant j3. Certain solutions for a particular 
set of initial parameters with a constant f3 do not exist when the variation of f3 
is incorporated appropriately. We model the viscous sub-Keplerian accretion disk 
with a nonzero component of advection and pressure gradient around black holes 
by preserving the conservations of mass, momentum, energy, supplemented by the 
evolution of j3. By solving the set of five coupled differential equations, we obtain 
the thermo-hydrodynamical properties of the flow. We show that during infall, (3 of 
the flow could vary upto ~ 300%, while 7 upto ~ 20%. This might have a signif- 
icant impact to the disk solutions in explaining observed data, e.g. super-luminal 
jets from disks, luminosity, and then extracting fundamental properties from them. 
Hence any conclusion based on constant 7 and (3 should be taken with caution and 
corrected. 
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1 Introduction 



Gas flows in the vicinity of a compact object, particularly black hole and 
neutron star, are expected to be highly relativistic. However, far away of it, 
when infalling matter comes off, e.g., a companion star, flows are rather non- 
relativistic. Therefore, during accretion of matter from a large distance, flows 
should transit from the non-relativistic to relativistic regime while reaching 
the vicinity of the black hole horizon or the neutron star surface, when flow is 
expected to be sub-Keplerian in nature. Indeed, it is generally believed that 
the velocity and temperature at the inner edge of accretion disk are very high. 

Not only the cases in accretion, highly relativistic flows are involved in astro- 
physical jets from galactic, extragalactic black holes, gamma-ray bursts etc. 
(e.g. Zensus 1997; Mirabel & Rodriguez 1999; Meszaros 2002), with the veloc- 
ity 0.9—0.98 times the speed of light. The collimated dipolar outflows emerging 
from deep inside collapsars, according to the collapsar model of gamma-ray 
burst (Woosley 1993), are expected to achieve a Lorentz factor more than 100. 
As the accretion and outflow/jet are expected to be coupled (influencing each 
other; Bhattacharya, Ghosh & Mukhopadhyay 2010), any change/transition 
(e.g. as mentioned above) in the accretion flow would influence the jet and 
hence the underlying inferences. 

When the velocity varies from the non-relativistic to relativistic (ultra- 
relativistic) regime, the corresponding polytropic index (7) in the equation 
of state (EOS), i.e. the ratio of gas to total pressure (0) in the present con- 
text, should not remain constant. Note that 7 and (3 can be approximated as 
constant only if the flow remains non-relativistic or ultra-relativistic through- 
out. Therefore, the existing solutions for such flows with a fixed 7 need to be 
verified and corrected if necessary with a relativistically correct EOS (see e.g. 
Chandrasekhar 1938). 

Most of the studies of accretion disks around compact objects have approx- 
imated 7 to be constant throughout the flow, which cannot be the correct 
description for the reasons described above. Blumenthal & Mathews (1976) 
introduced, for the first time to best of our knowledge, variable 7 and pro- 
posed an appropriate EOS in obtaining solutions for spherical accretion and 
wind around Schwarzschild black holes. Recently, Meliani et al. (2004) used 
their EOS and obtained the general solutions for the relativistic Parker winds 
incorporating varying 7 according to the nature of the flow. Ryu et al. (2006) 
proposed a new EOS suitable for non-relativistic and relativistic regimes with 
varying 7, to implement in the simulations of relativistic hydrodynamics. Very 
recently, Chattopadhyay & Ryu (2009) showed that composition of the fluid 
plays an important role to determine the value of 7 and then solutions of invis- 
cid spherical flows around black holes. This would have impact on a hot accre- 
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tion disk as well which was shown to exhibit significant nuclear burning, pro- 
ducing the variation of compositions in the disk (Chakrabarti & Mukhopad- 
hyay 1999, Mukhopadhyay & Chakrabarti 2000). Earlier, while Chen & Taam 
(1993) discussed structure and stability of the transonic optically thick accre- 
tion disks without restricting the gas and radiation content therein, did not 
show the actual variation of the gas/radiation content as a function of radial 
coordinate and then its impact to the solutions. 

The present work introduces a self-consistent variation of 7 in solving vis- 
cous flows in the accretion disk around black holes. We essentially concentrate 
on the region of the flows which is expected to be sub-Keplerian in nature. 
Therefore, following previous work (Muchotrzeb & Paczynski 1982; Abramow- 
icz et al. 1988; Chakrabarti 1996; Mukhopadhyay & Ghosh 2003; Rajesh & 
Mukhopadhyay 2010a) we model the flow from the Keplerian to sub-Keplerian 
transition radius (x t ) to the black hole event horizon (x + ). We plan to con- 
centrate, in particular, the class of flows which is not geometrically very thick, 
yet sub-Keplerian in nature including a nonzero component of advection and 
gradient of pressure in general. Hence, the underlying equations can be in- 
tegrated/averaged vertically without losing important physics, in obtaining 
the solutions. Therefore, following previous work (Muchotrzeb & Paczynski 
1982; Abramowicz et al. 1988; Chakrabarti 1996; Mukhopadhyay & Ghosh 
2003; Rajesh & Mukhopadhyay 2010a, 2010b), we consider vertically inte- 
grated infall consisting of gas and radiation. While the flow around Xt might 
be non-relativistic with a high 7 (depending on the size of the disk), as it 
advances towards x + , 7 must be varying. We plan to investigate, how the self- 
consistent variation of 7 over the disk radii affects the established solutions 
with a fixed 7. 

In the next section we establish the model equations describing flows. Subse- 
quently, we discuss solutions in §3 and finally summarize the results in §4. 



2 Model Equations 

The basic equations describing conservations of mass and momentum in the 
vertically integrated disk remain same as of earlier works which assumed 7 and 
then (3 to be constant throughout (e.g. Abramowicz et al. 1988; Chakrabarti 
1996; Rajesh & Mukhopadhyay 2010b). Abramowicz et al. (1988) also showed 
that the geometrically slimmer flow is stable for f3 > 0.4. For simplicity, fol- 
lowing e.g. Narayan & Yi (1994) and Chakrabarti (1996), we assume that the 
heat radiated out to be proportional to the heat generated by viscous dissipa- 
tion. Hence, following the same authors, we define a parameter / determining 
the fraction of matter advected in. Therefore, we recall the equation describ- 
ing conservation of energy in the vertically integrated disk from Chakrabarti 
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(1996) and Mukhopadhyay & Ghosh (2003). However, in reality, / should be 
determined, rather fixing, selfconsistently as done by Rajesh & Mukhopadhyay 
(2010a) very recently for hot optically thin flows around black holes. 

For the present purpose, all the above mentioned equations of conservation 
are supplemented by an equation describing the variation of /3 as a func- 
tion of radial coordinate, assuming a single temperature fluid. Earlier authors 
(e.g. Chen, Abramowicz & Lasota 1997) already formulated single tempera- 
ture optically thin disk flows including bremsstrahlung cooling for constant 
j3. Therefore, without repeating the description of equations of conservation, 
assuming the flow to be a mixture of perfect gas and radiation (Abramowicz 
et al. 1988; Chakrabarti 1996; Rajesh & Mukhopadhyay 2010a,b), we recall 
straight away 



_ -Pgas _ Pgas p _ pkT ^ 

P Prad ~\~ Pgas pi m p 

where k is the Boltzmann constant, m p the mass of a proton, pi the mean 
molecular weight of ions, and P, p, T are respectively the total pressure, 
density, temperature in K of the flow. Now for an optically thick flow, as was 
initiated in the limit of geometrically slim disk by Muchotrzeb & Paczynski 
(1982), Abramowicz et al. (1988), 



-Prad = -Pbb = pressure from blackbody radiation = — — , (2) 

where a is the radiation constant. On the other hand, for the optically thin 
hot flows (e.g. Narayan & Yi 1995; Mandal & Chakrabarti 2005; Rajesh & 
Mukhopadhyay 2010a,b) in the presence of bremsstrahlung as radiation mech- 
anism, which are presumably geometrically thicker than that of Muchotrzeb 
& Paczynski (1982), Abramowicz et al. (1988), 



-Prad = -Pbr = pressure from bremsstrahlung radiation 

~2x 10" 18 p 2 h T 1/2 (1 + 4.4 x 10~ 10 T) M (3) 

in CGS unit, where h is the half-thickness of the flow in units of GM/c 2 , G the 
Newton's gravitation constant, M the mass of the black hole, c the speed of 
light. Note that, in the first approximation (Rajesh & Mukhopadhyay 2010b), 
we do not consider other radiation mechanisms which could be effective in 
optically thin flows. In a future work, such effects e.g. synchrotron radiation, 
inverse-Comptonization could be included, once the present results confirm 
such a study worth pursuing. With the choice of 7P = pc 2 (Mukhopadhyay & 
Ghosh 2003), we further write 
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, 8-3/3 kT 



s 6-3/3/im p /3' 1 ; 

where c s is the sound speed of the flow. 

2.1 Optically thick flows 



First we set up the model equations for optical thick and geometrically slim 
flows. Eliminating T from eqns. ([I]), (J2J) and (j!J) we obtain 



jz jg + 6(l-^) ^, +C g_ 1) ljg = (5) 
ax c s ax pax 



where 



«4+3(l-ffl(l- (8 _ 3/j) 6 (6 _ 3/j) ) (6) 

and x is the radial coordinate in units of GM/c 2 . As /3 and 7 are related by 
(Narayan & Yi 1995; Ghosh & Mukhopadhyay 2009) 



from eqn. ([6]) one can easily obtain d^/dx as well. Now replacing P, using 
eqns. (CQ) and (jlj), from the momentum and energy conservation equations and 
then eliminating p using the mass conservation equation (see Mukhopadhyay 
& Ghosh 2003, for details), we obtain 



7 1\ d$ | 1 / L Wc, = 3 | 7 A 2 p \ J_dF 



Kc 2 "d I dx Co \K ) dx 2x Kc 2 \ x 3 I 2F dx 



A— + dB— + ^ InZ dX - 2 f InZX 3M + M dF / 9 n 
dx dx x dx x 2 2 x 2F dx ' 



.7; 

7) 



KZ-(1 + K)ad 2 } ^---(L-K + 1)(Z- ad 2 )^- + ^ 



dx c s dx dx 

\-\ K )z +{ K-^ + m ^^), (10) 
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where i? and A are the radial velocity and specific angular momentum of the 
disk respectively, / is kept constant for a particular flow (Mukhopadhyay & 
Ghosh 2003; Rajesh & Mukhopadhyay 2010b); / — > 1 corresponds to the 
advection of entire energy in and / — > to no advection, a the viscosity 
parameter, F the pseudo- Newtonian gravitational force due to a rotating black 
hole (Mukhopadhyay 2002), the functions A, B, Z, K, L are given by 



K 



7 r 3 - 1 ' 7 r 3 - 1 

,, 6Q8-1) 

(6 - 3/3)(8 - 3/3)i?' 



2 + 



36(1-/3) 
#(6-3/3X8-3/3)' 



(11) 



where 7 n , I n +i, Ti, r3 have their usual meaning, discussed by the previous au- 
thors in detail (e.g. Mukhopadhyay & Ghosh 2003; Rajesh & Mukhopadhyay 
2010a,b). 

Now combining eqns. ©, ©, fflQl) we obtain 



— = — , where 



(12) 



/i 



" 3 






1 dF~ 




2x " 






~ 2F~dx 





'Z - a 2 tf 



2m2\ 



c, Vtf 



2fI n ZX 3 Ad AddF\ 
x~ 2 2~ 2F dx ) 
I n+1 cl 1 dF\ 



- 2 K)Z + (K-l)i-a^ + xa ^ 



(13) 



1 

c.. 



(L-K + l)(Z-ad 2 



K 



-1)1 Ad - 



KZ-(l + K)ot& 
d 



(14) 



We strictly follow the methodology adopted by Rajesh & Mukhopadhyay 
(2010a) in solving eqn. f[T2~j) . Then following Mukhopadhyay & Ghosh (2003) 
and Rajesh & Mukhopadhyay (2010a,b) the outer boundary x t , where the flow 
becomes purely sub-Keplerian, is specified with the condition X/\k = 1 (Xk 
being value of Keplerian A). As matter advances from xt to x + , it reaches the 
critical radius x = x c when /i = f 2 = 0. Hence from / 2 = 0, the algebraic 
equation of Mach number (M ac ) at x c is given by 
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AM* c + BM 2 ac + C = 0, where 



(15) 



L-K + V x L-K + l 
K(T 3 -1) + K 




( 



| - 1) 




# 1 In 



l)(T 1 -K) + (L-K + T l ) 



(L 



K + l)a 2 fl n+l 





Then M ac can be easily computed from eqn. (TTBT) for given a and /, once 7 
and therefore (3 is known at x c , say 7 C and therefore (3 C . Subsequently, c s at 
x c , say c sc , can be obtained from fi = 0, once x c and the specific angular 
momentum of the disk at this radius A c are obtained. 

Now using the conservation of mass, eliminating dp/dx from eqn. ([5]), we 
obtain 



Hence, by solving eqns. (|T2l) and (jTTj) . along with any two of eqns. ([8]), ([9]), 
ffTUl) simultaneously, we obtain c s , A, (3 (and then 7) as functions of x from 

a; t to x + . 

2.2 Optically thin flows 

In presence of bremsstrahlung radiation for optically thin flows, from eqns. 
(CED, (ED and gl) we obtain 



dp _ 1 [ 7(1 - /3) cfc s (1 - /3) eft? 



<ic i? c s dx $ dx 



+ (1 





(18) 



where 



R 



(3 2 {(3 (8 — 3/3)(6 - 3/3) 



11-/3/1 6 




7 



The remaining equations are same as in an optically thick flow, with R given 
by eqn. (fT9~|) . Subsequently, as in §2.1, eliminating dp/dx from eqn. (fTBj) . we 
can obtain d(3/dx in terms of dd/dx and dc s /dx and hence solve for c s , A, /3 
(and then 7). 



3 Solution 

First we discuss the solution for the optically thick flows which are geomet- 
rically slimmer. Subsequently, we address the solution for the optically thin 
and geometrically thicker flows. 

3. 1 Optically thick flows 

We consider typical cases of accretion flow with (3 at x t being 0</3 t <l; extreme 
radiation dominated to extreme gas dominated flows at x t . Note that the 
general criterion for thermal stability of the flow is /3 > 0.4 (Abramowicz 
et al. 1988). We also assume that a part (/) of the energy dissipated in the 
infalling matter advects in{H However, note importantly that the value of / 
is not arbitrary, rather should be chosen according to the nature of the flow; 
whether it is gas dominated or radiation dominated, optically thin or thick. 
The flow consisting of significant radiation component should be radiatively 
more efficient and hence should have a small / compared to a flow with a lesser 
content of radiation. However, the flow with significant radiation around a 
rotating black hole might have a larger / compared to its non-rotating counter 
part, as the radial velocity of infall matter is higher for a rotating black hole 
with a shorter infall time scale. 

Figures [T^i,b show how 7 and /3 decrease as the flow advances from xt to 
x + . When p t ~ 1, the flow is extremely gas dominated at Xt and 74 ~ 5/3. 
However, as the flow advances in, the gas flow starts becoming relativistic, 
rendering the increase of apparent content of radiation which decreases (3. At 
x + , f3 decreases by ~ 30%, while 7 by ~ 9%. However, as (3 t and then j t 
decreases, their variation decreases as well. For fl t ~ 0.33, its variation during 
infall is only ~ 7%, while for a smaller fit variation is negligible. Therefore, for a 
gas dominated flow, the solutions of accretion flows, particularly for vertically 
integrated geometrically thin/slim disks, with a constant /3 and then 7 as 
discussed in the existing literature, appear to be incorrect which need to be 
revised. 

2 Qzvd = Q + ~ Q = f Q + 1 when Q av( j, Q + and Q~ respectively the amount of 
energy advected in, dissipated and radiated out. 



8 



Figures [Tfc,d show how the variation of affects the disk thermo- 
hydro dynamics, particularly for gas dominated flows. In general, as decreases 
at a large distance from the black hole, the content of radiation in the flow 
increases which leads to a smaller $ therein. This is because, the increase of 
radiation makes the flow to be radiatively more efficient, rendering the disk 
to be geometrically thinner and then more centrifugally dominated. Hence, 
ignoring the variation of 0, when t and C are high, would make the disk to 
be gas pressure dominated geometrically thick and quasi-spherical through- 
out which could reflect a wrong picture. However, at the inner region of the 
disk, d profiles for all merge. This is because the flow therein is practically 
controlled by the black hole's gravity. Interestingly, while the flow with a con- 
stant (= 0.75) would have a physical solution (dot-dashed line) for A c = 3.1, 
the variation of (with C = 0.75), when other initial parameters remain un- 
changed, leads to an unphysical solution (solid line) with an O-type critical 
point which does not continue from x t to x + . However, if A c decreases slightly 
(to , e.g., 3.05), then the flow with varying again exhibits a physical solution 
(long-short-dashed line). Figure [Tfe shows how the relation between P, p and 
c s , i.e. EOS, changes with the variation of in the flow. Clearly dP/dp ~ Pj ' p 
(which is practically true for the present purposes) does not vary linearly with 
c^, unlike the cases modeled in the existing literature based on constant 7. 

In Fig. [2] we compare the change of variation of 0, 7 and then corresponding 
with the change of a. It is found that the effect of variation of on disk 
thermo-hydrodynamics is stronger for a lower viscosity. If C is fixed, then is 
higher and its variation is steeper for a lower a at a large distance. Higher the 
0, higher the 7 in the flow is, inducing a quasi-spherical gas flow with a faster 
radial infall which may lead to an unstable/unphysical flow solution at a low 
a. In other words, a lower a corresponds to a lower rate of energy-momentum 
transfer in the disk. Hence, to describe a steady infall, has to be larger, 
particularly at a large x, which implies a quasi-spherical faster inflow of gas. 
Therefore, the assumption of a constant and then 7 (~ 1.533) would reflect 
an incorrect solution of the disk flow, particularly for a gas dominated flow. 
Figures[2b,d show that physically interesting solutions, exhibiting flows coming 
from x t to x + , require a smaller A c than that considered in the literature 
based on a constant 0. This supports quasi-spherical flows. Figure[2fe compares 
solutions with different sets of C and a, and shows that the physical solution 
is achieved for a smaller C , when a is lower. 

From Fig. [3] we understand how the variation of Kerr parameter affects the 
solutions with varying 0. First we consider the maximum possible A c obtained 
with constant (Rajesh & Mukhopadhyay 2010a) for the respective cases of 
a. Interestingly, while the physical solution around a co-rotating (prograde) 
black hole is available (solid line) with varying for the same set of parameters 
as of constant 0, A c has to be smaller for a counter-rotating (retrograde) 
black hole (dotted and dashed lines). This is because a high t renders a 
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quasi-spherical gas flow, which requires a smaller A allowing the matter to fall 
inwards. Note that while the maximum possible A c for a co-rotating black hole 
itself is smaller, for a counter-rotating black hole it is much larger. Hence, the 
allowed range of A c is shrunk for a counter-rotating black hole compared to 
that inferred from the analysis based on a constant /3. For clarity, in Figs. |3b,d 
we also show the $ and A profiles for a = 0.998 and —0.998 (respectively dot- 
long-dashed and dashed-long-dashed lines) with a constant for maximum 
possible values of A c . Figured shows how the EOS changes with the variation 
of a, when (3 varies in the flow. Clearly P/p does not vary linearly with c 2 s , 
unlike the cases with constant 7. Note that dashed and dashed-long-dashed 
lines practically overlap each other, implying that the solution with a constant 
(3 is recovered by decreasing A c in the case of varying f3. 

We have already seen in Fig. [I] that the variation of /3 is negligible for the 
radiation dominated flows. This is because a small corresponds to a radia- 
tively efficient geometrically thinner flow, which effectively corresponds to a 
relativistic flow of radiation throughout. Therefore, as the matter advances, 
the black hole's gravity has nothing to affect especially. However, close to the 
black hole, 1? increases noticeably and the corresponding infall time scale of the 
matter decreases significantly, rendering a quasi-spherical flow which increases 
(3 by 4%, as shown in Fig. 0} 

Figure shows that the variation of / keeping other parameters intact does 
not affect hydrodynamic properties, e.g. Mach number, significantly. However, 
as / increases, the flow tends to become radiatively inefficient, rendering an 
increase in (3 significantly close to the black hole. Therefore, any (observed) 
property related to the inner accretion flow must be influenced by / and hence 
the variation of (3. 



3.2 Optically thin flows 

There is a stronger variation of (3 and then 7 in optically thin flows, compared 
to its optically thick counter parts. As seen in Fig. El while for a high (3 t , the 
variation of (3 could be ~ 100% and of 7 be ~ 15%, for a low (3 t , the respec- 
tive variations are much higher, ~ 300% and ~ 20%. In all the cases depicted 
here, as A decreases around x t (see Fig. [6H), rendering the flow to be quasi- 
spherical, (3 increases as the flow advances in. Subsequently, centrifugal force 
starts increasing, which renders the increase of A and rate of cooling, which 
decreases (3 and 7. However, close to the black hole, due to the dominance 
of gravitational force, matter falls very fast, rendering the flows, irrespective 
of /3t, to be quasi-spherical with high j3. Therefore, unlike our conventional 
thought that highly relativistic flow in the vicinity of a black hole may exhibit 
low (3 and then 7 as seen in optically thick flows, in optically thin cases with 
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bremsstrahlung radiation the situation is different. This is because the ineffi- 
cient/incomplete cooling process (Rajesh & Mukhopadhyay 2010a) keeps the 
flow to be very hot and then quasi-spherical until inner edge when the infall 
time is very short. It is very clear from EOSs shown in Fig. that P/p does 
not vary linearly with c 2 s , unlike the conventional cases. 

Figure [7] shows that the self-consistent inclusion of variation of ft and then 
7, for both the gas and radiation dominated flows, renders the unphysical 
solutions (solid lines), when the corresponding solutions with constant ft and 
then 7 (dotted lines) appear physical. It is very clear that at around x = 50, the 
value of 7 varies sharply with radius, rendering the flow to be gas dominated 
quasi-spherical with an unphysical O-type critical radius. However, the flow 
with a smaller A c (dashed lines), keeping other parameters unchanged, exhibits 
the physical solution, even with varying ft. This is because the decrease of A 
favors the quasi-spherical nature of gas dominated flow at around x>50, giving 
rise to a physical solution. In general, the specific angular momentum of the 
flow, more specifically A c , is smaller in a realistic flow with varying ft and 7 
than that with the idealistic choice of constant ft and 7. 

Figure [8] shows that the variation of / keeping other parameters intact neither 
affects hydrodynamic properties nor ft profiles. However, at a very high / (— > 1, 
advection dominated flow), other parameters should not be kept same and ft 
must be larger. Such a flow is represented by long-dashed lines. However, this 
still does not alter the Mach number profile significantly. Note also that the 
flow with ft = f = 1 strictly at a large distance does not exhibit any variation 
of ft even at a smaller radius, which renders the corresponding 7 = 5/3, to be 
ratio of specific heats throughout. Therefore a strict ADAF solution remains 
unaltered. However, Narayan Sz Yi (1994) in their ADAF chose 7 = 1.5, while 
7 = 5/3 corresponds to a Bondi flow. Moreover, for all the practical flows 
ft, f < 1 and hence a variation of ft should be accounted in the model, in 
particular to explain luminosity, as the present work intends to do. 



4 Summary 

We have investigated the sub-Keplerian general advective accretion flows 
(GAAF; Rajesh & Mukhopadhyay 2010a) around black holes allowing evo- 
lution of the gas and radiation content into flows self-consistently. We have, 
in one hand, considered radiation to be optically thick blackbody in the geo- 
metrically thin/slim vertically integrated disks. On the other hand, we have 
also considered optically thin geometrically thicker flows in the assumption of 
bremsstrahlung radiation. Hence, as matter advances towards a black hole, the 
ratio of the gas pressure to total pressure and the corresponding polytropic 
index vary significantly with the radial coordinate. We have found that in 



11 



Table 1 

Parameter sets used for optically thick solutions: M = 10, m is the accretion rate 
in Eddington units 



Br 


IC 




f 




CL 




777 


0.75 


1.533 


3.1 


0.5 


5.8 


o 


0.01 


0.01 


0.75 


1.533 


3.1 


0.5 


5.8 


o 


0.0001 


0.01 


0.75 


1.533 


3.05 


0.5 


5.8 


o 


0.01 


0.01 


0.75 


1.533 


3.05 


0.5 


5.8 


o 


0.0001 


0.01 


0.745 


1.531 


3.1 


0.5 


5.8 





0.01 


0.01 


0.745 


1.531 


3.1 


0.5 


5.8 





0.0001 


0.01 


0.7 


1.513 


3.1 


0.5 


5.8 





0.0001 


0.01 


0.75 


1.533 


1.7 


0.5 


3.5 


0.998 


0.01 


0.01 


0.75 


1.533 


4 


0.5 


7.5 


-0.998 


0.01 


0.01 


0.75 


1.533 


3.95 


0.5 


7.5 


-0.998 


0.01 


0.01 


0.5 


I AAA 


3.22 


0.5 


5.8 





0.01 


0.01 


0.5 


I AAA 


3.22 


0.5 


5.8 





0.01 


0.01 


0.33 


1.399 


3.2 


0.1 


5.8 





0.01 


1 


0.33 


1.399 


3.2 


0.3 


5.8 





0.01 


1 


0.33 


1.399 


3.2 


0.7 


5.8 





0.01 


1 


0.33 


1.399 


3.2 


0.9 


5.8 





0.01 


1 


0.05 


1.342 


3.2 


0.3 


5.8 





0.01 


1 


0.05 


1.342 


1.8 


0.5 


3.5 


0.998 


0.01 


1 


constant 0/j 


0.75 


1.533 


3.1 


0.5 


5.8 





0.01 


0.01 


0.75 


1.533 


3.1 


0.5 


5.8 





0.0001 


0.01 


0.75 


1.533 


1.7 


0.5 


3.5 


0.998 


0.01 


0.01 


0.75 


1.533 


4 


0.5 


7.5 


-0.998 


0.01 


0.01 



several occasions, the accretion solutions as obtained with a constant (3 by the 
earlier models do not exist when f3 is allowed to vary self-consistently. It has 
a very important implication when most of the existing models of accretion 
disks assume (5 and 7 to be constant throughout the flow. For example, an ad- 
vection dominated class of solution, namely ADAF (Narayan & Yi 1994, 1995; 
Quataert & Narayan 1999), has been modeled based on gas dominated flows 
with a constant 7 (= 1.5) throughout, when 7 is chosen to be ratio of specific 
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Table 2 

Parameter sets used for optically thin solutions: x c = 5.8, a = 0, a = 0.01 (except 
for the case with f3 c = 0.95, when a = 0.001), M = 10, m is the accretion rate in 
Eddington units 





7c 


A c 


/ 


m 


0.75 


1.533 


3.14 


0.5 


0.01 


0.75 


1.533 


3.1 


0.5 


0.01 


0.75 


1.533 


3.1 


0.1 


0.01 


0.75 


1.533 


3.1 


0.8 


0.01 


0.95 


1.635 


2.8 


0.98 


0.01 


0.5 


1.444 


3.22 


0.5 


0.01 


0.5 


1.444 


3.2 


0.5 


0.01 


0.33 


1.399 


3.25 


0.5 


1 


0.33 


1.399 


3.2 


0.5 


1 


0.05 


1.342 


3.27 


0.5 


1 


0.05 


1.342 


3.25 


0.5 


1 


constant 


0.75 


1.533 


3.14 


0.5 


0.01 


0.5 


1.444 


3.22 


0.5 


0.01 


0.33 


1.399 


3.25 


0.5 


1 


0.05 


1.342 


3.27 


0.5 


1 



heats. Other models (e.g. Abramowicz et al. 1995; Chakrabarti 1996; Rajesh 

6 Mukhopadhyay 2010a), describing gas dominated optically thin flows, as 
well have assumed 7, which is polytropic index therein, to be constant. In 
the present work, the impact of varying 7 has been particularly demonstrated 
in the framework of later model. Note, however, that polytropic index turns 
out to be a ratio of specific heats only when the flow is of pure gas with 

7 = 5/3. The results argue for correction to the earlier models based on con- 
stant /3, when particularly (3 varies more than 100% during the infall of an 
optically thin gas dominated matter. Hence, the present results imply that the 
entire branch of optically thin hot solutions of accretion disk needs to be re- 
investigated including appropriate variation of 7 and (3. However, for simplicity 
we have considered the optically thin radiation arising due to bremsstrahlung 
mechanism only, neglecting synchrotron, inverse-Comptonization which are 
likely to operate in the inner hot regime of the flow. However, inclusion of 
such effects would only favor the present results. In future, we have plan to 
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2 

log(x) 



0.02 0.04 

n 2 



0.06 



Fig. 1. Variation of (a) polytropic index, (b) ratio of gas pressure to total pressure, 
(c) Mach number, (d) ratio of specific angular momentum to the corresponding Ke- 
plerian value. In (a), (b), (c), (d) the solid, dotted, dashed, dot-long-dashed lines 
corresponding to the different j3 at transition radius are for j3 c = 0.75, 0.5, 0.33, 0.05 
respectively. In (c) and (d) long-short-dashed and dot-dashed lines correspond to 
j3 c = 0.75, but respectively with a smaller A c (= 3.05) than that of solid line 
(A c = 3.1) and with a constant 7 with A c = 3.1. (e) The density-pressure-sound 
speed relation (EOS), when the profiles shown by solid and dotted lines are the 
corresponding EOS of the solutions shown in (c) and (d) by solid and dotted lines 
respectively, while dashed line corresponds to a case with smaller A c than that of 
solid line. Other parameters are a = 0.01, / = 0.5, x c = 5.8, a = 0, M = 10. See 
Table 1 for detail. 



include such physics in studying very hot flows. This will have immediate in- 
fluence to the computed spectra (e.g. Yuan et al. 2003; Mandal & Chakrabarti 
2005) and subsequent inferences about sources. 
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Fig. 2. Variation of (a) polytropic index, (b) ratio of gas pressure to total pressure, 
(c) Mach number, (d) ratio of specific angular momentum to the corresponding 
Keplerian value, as functions of radial coordinate. In (a), (b), (c), (d) the solid and 
dotted lines are for a = 0.01 and 0.0001 respectively when A c = 3.1, and dashed 
and dot-long-dashed lines for a = 0.01 and 0.0001 respectively when A c = 3.05. In 
(c) and (d) long-short-dashed and dot-dashed lines are corresponding solutions with 
constant 7 (~ 1.533) for the cases of a = 0.01 and a = 0.0001 respectively when 
A c = 3.1. (e) Mach number as a function of radial coordinate, when solid and dotted 
lines are for a = 0.01 and 0.0001 respectively where (3 C = 0.745, and dashed line 
is for a = 0.0001 where j3 c = 0.7. Other parameters are / = 0.5, x c = 5.8, a = 0, 
M = 10. See Table 1 for detail. 
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Fig. 3. Variation of (a) polytropic index, (b) ratio of gas pressure to total pressure, 
(c) Mach number, (d) ratio of specific angular momentum to the corresponding 
Keplerian value. In (a), (b) the solid, dotted, dashed, dot-long-dashed lines respec- 
tively correspond to parameter sets a = 0.998, A c = 1.7; a = —0.998, A c = 4; 
a = -0.998, A c = 3.95; a = 0, A c = 3.05. In (c), (d) the solid, dotted, dashed, dot- 
long-dashed, dashed-long-dashed lines respectively correspond to parameter sets 
a = 0.998, A c = 1.7; a = -0.998, A c = 4; a = -0.998, A c = 3.95; a = 0.998, A c = 1.7, 
with constant f3 = 0.75; a = —0.998, A c = 4, with constant /3 = 0.75. The profiles 
shown by solid, dotted, dashed, dashed-long-dashed lines in (e) are the correspond- 
ing density-pressure-sound speed relation (EOS) of the solutions shown by the same 
lines in (c) and (d). Other parameters are a = 0.01, / = 0.5, M = 10. See Table 1 
for detail. 
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